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Abstract 

We study the impact of semi-annihilations XiXj ^ XkX, where Xi is any dark matter and X is 
any standard model particle, on dark matter phenomenology. We formulate minimal scalar dark 
matter models with an extra doublet and a complex singlet that predict non-trivial dark matter 
phenomenology with semi-annihilation processes for different discrete Abelian symmetries Zj^ , N > 
2. We implement two such example models with and Z4 symmetry in micrOMEGAs and work 
out their phenomenology. We show that both semi-annihilations and annihilations involving only 
particles from two different dark matter sectors significantly modify the dark matter relic abundance 
in this type of models. We also study the possibility of dark matter direct detection in XENONIOO 
in those models. 

1 INTRODUCTION 

The origin of dark matter of the Universe is not known. In popular models with new particles beyond 
the standard model particle content, such as the minimal supersymmetric standard model, an addi- 
tional discrete Z2 symmetry is introduced |1]. As a result, the lightest new Z2-odd particle, x, is stable 
and is a good candidate for dark matter. The phenomenology of this type of models has been studied 
extensively. 

The discrete symmetry that stabilises dark matter must be the discrete remnant of a broken gauge 
group [2], because global discrete symmetries are broken by gravity. The most natural way for the 
discrete symmetry to arise is from the breaking of a U{l)x embedded in a larger gauge group, e.g. 
50(10) [3]. The latter contains gauged i? — L as a part of the symmetry, and the existence of dark 
matter can be related to the neutrino masses, leptogenesis and, in a broader context, to the existence 
of leptonic and baryonic matter [4||6] . 

Obviously, the discrete remnant of U{l)x need not to be Z2 - in general it can be any Z^v Abelian 
symmetry. The possibility that dark matter may exist due to Zjv, > 2, is a known [7]-|15|, but 
much less studied scenario Model independently, it has been pointed out in Ref. [is] that in Zjy 
models the dark matter annihilation processes contain new topologies with different number of dark 
matter particles in the initial and final states - called semi- annihilations -, for example xx o x*X, 
where X can be any standard model particle. It has been argued that those processes may significantly 
change the predictions for the dark matter relic abundance in thermal freeze-out. Furthermore, an 



^Phenomenology of Za-symmetric dark matter in supersymmetric models has been studied in Refs. 7][To and in extra 
dimensional models in Refs. i8||9|. 



enlarged discrete symmetry group makes it possible to have more than one dark matter candidate. In 
this case, annihilation processes involving only particles from the dark sectors, leading to the assisted 



freeze-out mechanism, can also influence the relic abundance of both dark matter candidates 16, 17 



The assisted freeze-out mechanism in the case of a Z2 x Z2 symmetry was discussed in 17 . However, 
no detailed studies have been performed that compare dark matter phenomenology of different 
models. This is difficult also because presently the publicly available tools for computing dark matter 
relic abundance do not include the possibility of imposing a Zfq discrete symmetry instead of a Z2. 

The aim of this work is to formulate the minimal scalar dark matter model that predicts different 
non-trivial scalar potentials for different Zfq symmetries and to study their phenomenology. In par- 
ticular we are interested in quantifying the possible effects of semi-annihilation processes xx -f-)- x*X 
as well as of annihilation processes involving particles from two different dark sectors on generating 
the dark matter relic abundance. In order to perform quantitatively precise analyses we implement 
minimal and Z4 symmetric scalar dark matter models that contain one singlet and one extra dou- 



blet in micrOMEGAs 18, 19 . Using this tool we show that, indeed, the semi-annihilations and the 
annihilations between two dark sectors affect the dark matter phenomenology and should be taken into 
account in a quantitatively precise way in studies of any particular model. 



2 LAGRANGIANS 

2.1 Zjq symmetry 

Under an Abelian Zj<; symmetry, where is a positive integer, addition of charges is modulo A^. Thus 
the possible values of Z^ charges can be taken to be 0, 1, ... ,A^ — 1 without loss of generality. A 
field with Zj^ charge X transforms under a Zj^ transformation as </> — )• (\)^ where = 1, that is 
Lo = exp(i27r/A'). 

A Zjv symmetry can arise as a discrete gauge symmetry from breaking a U{l)x gauge group 
with a scalar, whose X-charge is N [2|[4]. For larger values of A^, the conditions the Zn symmetry 
imposes on the Lagrangian approximate the original U{1) symmetry for two reasons. First, assuming 
renormalisability, the number of possible Lagrangian terms is limited and will be exhausted for some 
small finite A^, though they may come up in different combinations for different values of A^. Second, 
if the symmetry arises from some U{l)x, the X-charges of particles cannot be arbitrarily large, 
because that would make the model nonperturbative. If A^ is larger than the largest charge in the 
model, the restrictions on the Lagrangian are the same as in the unbroken U{1). 

We shall see below that in spite of the large number of possible assignments of Z]\f charges to the 
fields, the number of possible distinct potentials is much smaller. 



2.2 Field content of the minimal model 

In order to study the impact of different discrete Z^ symmetries on dark matter phenomenology, 
the example model must contain more than one neutral particle in the dark sector. The minimal 
dark matter model with such properties contains, in addition to the standard model fermions and the 
standard model Higgs boson Hi, one extra scalar doublet H2 and one extra complex scalar singlet S [5]. 
In the case of Z2 symmetry, as proposed in [s], those new fields can be identified with the well known 



inert doublet H2 |20]^23j and the complex singlet S [24-28 . The phenomenology of those models is well 
studied. However, when both the doublet and singlet are taken into account, qualitatively new features 
concerning dark matter phenomenology, electroweak symmetry breaking and collider phenomenology 
occur [5|[6| [29] - [3T] . The field content of the minimal scalar Z^r model is summarised in Table [Tj 
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Table 1: Scalar field content of the low energy theory with the components of the standard model 
Higgs Hi in the Feynman gauge. The value of the Higgs VEV is v = 246 GeV. 

Field 5C/(3) SU{2)l Y/2 Q = T^ + Y/2 



H 



1 - I^ D+fe+iG" j ^ ^ 1 1 ) 2 



Ho — rrO , • 1 2 



S = ^H+iSA I I 

v2 




1 



1 



2.3 Constraints on charge assignments 

The assignments of Zj^ charges have to satisfy 

Xs>0, 

-Xi + Xl + = mod N, (1) 
-Xq + Xi+Xd = Q mod N, 
-Xg -Xi+Xu = mod N. 

The first and second conditions arise from avoiding the \Hi\'^S term and Yukawa terms for H2, respec- 
tively, and the rest from requiring Yukawa interactions between Hi and standard model fermions. The 
choice of charges for standard model fermions, the standard model Higgs Hi, the inert doublet H2 
and the complex singlet S must be such that there are no Yukawa terms for H2 and no mixing between 
Hi and H2: only annihilation and semi-annihilation terms for H2 and S are allowed. While we will 
see below that there are many assignments that satisfy Eq. ([T]), in each case it was possible to find an 
assignment with the charges of standard model fields set to zero: -'^g/^M^d^e,! = 0. 

All possible scalar potentials contain a common piece because the terms where each field is in pair 
with its Hermitian conjugate are allowed under any Zjy and charge assignment. We denote it by Vc 
(the 'c' stands for 'common'): 

(2 \ 2 
l^il'-yj + film^ + X2\H2\^ + i^l\S\^ + Xs\S\^ ^2) 

+ \si\S\^\Hi\^ + Xs2\S\^\H2\^ + X3\Hi\^\H2\^ + Xa{hIH2){hIHi). 

2.4 The Z2 scalar potential 

There are 256 ways to assign the possible Z2 charges 0, 1 to the standard model and dark sector fields. 
Of these, 8 satisfy Eq. Q; among them, there are 2 different assignments to the dark sector fields: 
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Xs = Xi = 1, X2 = and Xi = 0, X2 = Xs = 1. Both give rise to the unique scalar potential 
V = V, + ^{S^ + 5t2) + ^ [{H\H,f + [HlH.f 

2.5 Z3 scalar potentials and particle content 

There are 6561 ways to assign 0, 1,2 to the fields. Of these, 108 satisfy Eq. ([T]); among them, there 
are 12 different assignments to the dark sector fields, giving rise to 2 different scalar potentials. The 
example potential we choose to work with (given by e.g. Xi = 0,X2 = Xs = 1) is 

I 2. (4) 

+ ^^{ShIHi + S^h\H2)^ 

which induces the semi-annihilation processes we are interested in. The second one is obtained from 
Eq. Q by changing 5-^5^ (with ^SH fJ-'sH ^^"^ -^512 A52i)- 

The following conditions are sufficient to have the global minimum of potential at electroweak 
vacuum with (5) = 0, (-^2) = 0: 

Ai,A2, A5,A5i, A52 > 0, (5) 

A3 + A4 > 0, (6) 

4A51A52 > A|i2 ) (7) 

As A3 + A4 
We use these conditions for our benchmark points. 

The last term in Eq. (|4]) induces a mixing between the down component of H2 and S. In terms of 
the mass eigenstates xi, X2, we have 

/ -iH+ \ 

H2 = \ , 5 = xi COS0 — a;2sin0. (9) 

\ xi sin + X2 cos J 

The dark sector of this model consists of 3 complex particles xi, X2, and with the Z3 charge of 
1. Taking the masses of xi, X2 and the mixing angle 9 as free parameters of the model, we get the 
following relations 

4 = Mlsm^e + Mlcos'e-Xsi'^, (10) 



, / , ,9 , ^9 N COS a sm (j ,^ , , 

^,SH = -MM^-Ml) — ^ — , (11) 

2 



V2v 

V 

~2 

The Ai and the mass of H'^ can be presented by formulas 



^2 = _(A4 + A3)— + M2^sin20 + M2^cos20. (12) 



- (13) 



Mh+ = \/^i + A3y. (14) 
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where M/^ is mass of SM Higgs. 



2.6 Z4 scalar potentials and particle content 

There are 65536 ways to assign 0, 1, 2, 3 to the fields. Of these, 576 satisfy Eq. Q; among them, there 
are 36 different assignments to the dark sector fields, giving rise to 5 different scalar potentials. Among 
those the only potential that contains semi-annihilation terms is 

Vk = V. + ^{S' + 5t4) + ^ UHlH.f + (HlH.f 

2 2 L J (15) 

+ ^{s^h\h2 + s^'Hlm) + ^(s'Hlm + s^'hIh,), 

invariant under e.g. the assignment of Z4 charges Xi = 0,X2 = 2,Xs = 1. 

The following conditions are sufficient to have global minimum of potential at electroweak vacuum 
with (5) = 0, {H2) = 0: 

Ai,A2,A5i,A52 > 0, (16) 
As-lA'sl > 0, (17) 
A3 + A4-IA5I > 0, (18) 

(IA512I + |A52l|)^ < AsiA52. (19) 

Our benchmark points considered below satisfy these conditions. 

The other four scalar potentials can formally be obtained from the Z2-invariant potential Eq. ([3]) 
by setting all the new terms added to Vc to zero, with the exception of the 1) A'^, ^j-sh, 2) A'^., fJ.'sjj, 3) 
l^'s^ A'^, A'^, X'g^, A'5.2, 4) /i'^,, A'^, A'^, A'^.^^, A'5.2, fJ-SH, fJ-'sH terms. 



The A5 term in potential (15) splits the down component of H2 into two real scalar fields with 
different masses, 

H2=(^ ^ . (20) 

Note that the complex scalar S does not mix with H2 because these fields have different Zn charges. 
As a result this model contains two dark sectors, the first one with the complex scalar S (the Z4 charge 
is 1), the second one comprising the complex scalar and the real scalars and ( the Z4 charge 
is 2). Any of the neutral particles with a non-zero Z4 charge can be a dark matter candidate. We will 
consider the masses of the neutral scalar particles, Ms, Mjjo and M^o, as independent parameters, 
then 



4 = Mj-Asiy, (21) 



A. ^ 'Ji^. (22, 

yZ 

fll = M^o - (A3 + A4 + A5)y, (23) 



^ (24) 
A. . \§. (25) 



5 



3 RELIC DENSITY IN CASE OF THE SYMMETRY 



3.1 Evolution equations 

Consider the Zs-symmetric theory. The imposed Zj, symmetry imphes, as usual, just one dark matter 
candidate. This is because the Z3 charges 1 and —1 correspond to a particle and its anti-particle. The 
new feature is that processes of the type xx — >■ where X is any standard model particle, also 

contribute to dark matter annihilation. The equation for the number density reads 



— = - - (n2 - nn) - 3Hn, (26) 

dt ^ 

where we use n = neq, H is the Hubble rate, and angular brackets mean thermal averaging. We define 



1 1 xx^x*X 

a, = (w""*^^^) + ^(w""^"*^) and a = , (27) 



2' ' 2 a. 



V 



which means that < a < 1. Here and in the following we use the notation, a^^^^ ^ = {x)a^^^^ ^). 
In terms of the abundance, Y = n/ where s is the entropy density, we obtain 

^ = -sa, _ aYY - (1 - a)Y^) (28) 
or, using the entropy conservation condition ds/dt = —3Hs, 

.dY /^,2 ..^ 



3H— =ay(Y^ -aYY -{l-a)Y) . (29) 



where Y = Y^q is the equilibrium abundance. We use standard formulae for H[T) and s{T) 32 that 



allow to replace the entropy evolution with the temperature one. To solve this equation we follow the 



usual procedure 18 ,32 . Writing Y = Y + Ay we find the starting point for the numerical solution of 



this equation with the Runge-Kutta method using 

dY — 

3H— = a^YAY (2 - a) , (30) 
ds 

where Ay ^ Y. This is similar to the standard case except that Ay increases by a factor 1/(1 — a/2). 
Furthermore, when solving numerically the evolution equation, the decoupling condition y^ ^ y is 
modified to 

y2 >ayy + (i-a)yl (31) 

This implies that the freeze-out starts at an earlier time and lasts until a later time as compared 



with the standard case. This modified evolution equation is implemented in micrOMEGAs [19 33 
Although semi-annihilation processes can play a significant role in the computation of the relic density 
the solution for the abundance depends only weakly on the parameter a, typically only by a few percent. 
This means in particular that the standard freeze-out approximation works with a good precision. 



3.2 Numerical results with micrOMEGAs 

Using the scalar potential defined in Eq. Q we have implemented in micrOMEGAs the scalar model 
with a Z3 symmetry. The scalar sector contains an additional scalar doublet and one complex singlet. 
The neutral component of the doublet mixes with the singlet, the lightest component xi is therefore 
the dark matter candidate, while the heavy component X2 can decay into xih, where h is the standard 
model-like Higgs boson. Because h can decay into light particles, X2 is unstable even if the mass 
difference between xi and X2 is small. Note that the doublet component of DM has a vector interaction 
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with the Z . This interaction is determined by the SU{2) x U{1) gauge group and leads to a large direct 
detection signal in conflict with exclusion limits, for example from XENONIOO |34|. The only way to 
avoid this constraint is to consider a DM with a very small doublet component, namely we have to 
assume that the mixing angle 

d < 0.025. (32) 

In the limit of small mixing, annihilation processes such as xix\ — t- XX where X stands for W, Z, h, are 
dominated by the Asi|5p|-ffi|^ term. The semi-annihilation process xixi — )■ x^h is mainly determined 
by a product of fig and A51 arising from the terms //'^(S*^ + 5'"'"'^)/2 and A5i|5'p|/fip in Eq.[2]and Eq.jij 
To illustrate a scenario where semi-annihilation channels contribute significantly and which predicts 
reasonable values for the relic density and the direct detection rate, we choose a benchmark point with 
the following parameters 



A2 


0.1 


Xs 


0.2 


As'12 


0.1 




150 GeV 


A3 


0.1 


A51 


0.05 


Mh 


125 GeV 




400 GeV 


A4 


0.1 


A52 


0.1 




80 GeV 


sin^ 


0.025 



Table 2: Benchmark point for Z3. 

For this point, the relic density is VLh? = 0.105. The dominant contribution to {U,h?)~^ is from semi- 
annihilation (54% for xixi — ?• /ix^) while the annihilation channels xix\ — >• WW, ZZ, hh give a relative 
contribution of 22%, 13% and 10% respectively. Fig. [T] illustrates the dependence of the relic density on 
the DM mass as compared to the relic density when semi-annihilation is ignored, (r2/i^)ann- Here all 
other parameters are fixed to their benchmark values. When Mx^ = 110 GeV, semi-annihilation with a 
Higgs in the final state is kinematically forbidden at low velocities. If Mx-^ increases, semi-annihilation 
plays an important role and il/i^ decreases rapidly due to the contribution of the channel xixi — )• hx\. 
Note that {^h?)a,nn also decreases when M^^ is such that the channel xix\ — )• hh is allowed. When 
Mx-^ approaches Mx^l'i'^ ^h'^ falls again because the semi-annihilation channel is enhanced due to X2 
exchange near resonance. 

The spin independent (SI) scattering cross section on nuclei as a function of the DM mass is 
illustrated in Fig. [l] (right panel). Here we average over dark matter and anti-dark matter cross 
sections assuming that they have the same density. The main contribution comes from the Z-exchange 
diagram because there is a xix\Z coupling Furthermore, one can easily show that the scattering 
amplitudes are not the same for protons and neutrons, with fp = (4sin^ 9\y — l)fn = — 0.075/„. Since 
the current experimental bounds on cjfp are extracted from experimental results assuming that the 
couplings to protons (fp) and neutrons (fn) are equal and the same as the couplings of to protons 
(fp) and neutrons (/n), we define the normalised cross section on a point-like nucleus 



35 



SI _ ^ , MnMx, Y ( [Zfp + [A-Z)Uf izfp + {A-Z)hf \ 



vr \Mn + MxJ \ A2 
This quantity can directly be compared with the limit on a^. 



In the inert doublet model with a Z2 symmetry 20 22 , a A5 term splits the complex doublet into a scalar and a 



pseudoscalar, when the mass splitting is small such coupling leads to inelastic scattering. 
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Figure 1: (Left panel) Qh'^ as a function of the dark matter mass for the benchmark point with semi- 
annihilation (solid line), and without semi-annihilation (dashed). (Right panel) <7||xe (solid). The 
experimental limit from XENONIOO 34 is also displayed (dashed). 



4 RELIC DENSITY IN CASE OF THE SYMMETRY 
4.1 Evolution equations 

In the case of a symmetry all particles can be divided into 3 classe^ {0,1,2} according to the value 
of their Z4 charges modulo 4. We can choose SM particles to have X^m = 0. We will use the notation 
(jO-bcd £qj, ^YiQ thermally averaged cross section for reactions ab — )• cd where a, b,c,d = 0, 1, 2 represent 
any particle with given X-charge. Let M^^^ and be the masses of the lightest particles of classes 
1 and 2 respectively. The lightest particle of class 1 is always stable and therefore a DM candidate. 
The lightest particle of class 2 is stable and can be a second DM candidate if < 2Mx^. Note that 
if > '2Mx^ , then X2 will decay before the freeze-out of xi and the relic density can be computed 
following the standard procedure. 

The equations for the number density of particles 1 and 2 read 



'^"'1 _ ^1100 {2 -2\ „1120 [2 -2"'2\ _1122 [2 ^2"-! 



\ 1120 / 2 -2'"^ \ 1122 2 2 "'l \ ou foA\ 

)-(7^ [ni-n^—]-a^ I - ng^g - 3i^ni, (34) 



^"■2 _ _2200 /'^2 -2\ , -'^_1120/„2 ^2^2\ 1^1210 



'-2 



2200 / 2 -2\ I 1120 / 2 -2 '"^ \ 1210 r - \ 



.^2211|^2_^2!^)_3^^^^ (35) 



^1 



where we use nj to designate the equilibrium number density of particle Xj. In a"^'^'^ all annihilation 
and coannihilation processes are taken into account. Here the semi-annihilation processes include all 
those, where 2 DM particles annihilate into one DM and one standard particle, specifically al^"^^ and 
o'l'^^^. These two cross sections are also described by the same matrix element. However, there is no 
simple relation between these two cross sections because one process is in the s-channel and the other 



^We take into account that 3 = — 1 mod 4, so the particle with X-charge 3 is the antiparticle of a particle with 
X-charge 1. 



8 



in the t-channel. In terms of the abundance, Yi = Ui/s, 

Q u '^^l ^1100 /ta2 a72\ , ^1120 -i^ ^1 \ , ^1122 /'ta2 v2^i\ /oa\ 

3-H = (^v [^1 -^i) +(^v I ^1 - ^2^ I + 0-^ I - >^ =2 I , (36) 
3H^ = °° {Yi - Yl) - lal'^' (y,' " >^2|^) + 1^1^''^ {Y, - ¥,) 

-'rf^'lYi-Yi^]. (37) 



Y 



Solving these equations we use standard formulas for entropy s{T) and the Hubble rate H{T) temper- 
ature dependence [32] that allow to replace the dependence on entropy with one on temperature. The 
thermally averaged cross section involving particles of different sectors can be expressed as 

I |-#ab->cd(\/i,cose)|2dcose, (38) 

pol. 

Yi{T) = ^^5,mfK2(^), (39) 

iei 

where ^ab^cd is the matrix element for the 2—7-2 process and Ki , K2 are modified Bessel functions of 
the second kind. For reactions which are kinematically open at zero relative velocity, depends slowly 
on temperature. Otherwise there is a strong exp(— AM/T) temperature dependence, where AM is the 



difference between the sums of the masses of outgoing and incoming particles. Equation (38) leads to 
relations between different cross sections 



YjYjai'^'^ = YKYLa^^''. (40) 

In particular it implies that, cr^^^^ = a^^'^^Y^ /Y2, where the abundance of incoming SM particles 
Yo = l. 



Introducing AYi = Yi — Yi, Eqs. (36) and (37) take a simple form 



AY 

= -C, + Aij{T)AYj + Qi,k{T)AYjAYk, (41) 



where 



a = 3H^, (42) 
as 

^ ^ , 2(^1100 + aii22 + aii20)Fi -(^^120 + 2^1122)2 

' -^1120^1 - 2all22y^ 2(cTf00+a2211)F2 + 0.5(^1210 + ^1120 yi)y^ I' ^ ^ 



2 



/ 1100+1122+1120 \ 

Qi = { \ 2211 1 ' (44) 



-a. 



-.1120 _ ^1122 lrrl210 



^2=1 ' . ^ .2200: .211 • (45) 



CT^'^ + a. 



(46) 
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At large temperatures we expect the densities of both DM components to be close to their equi- 
librium values. In general in micrOMEGAs [36] the equation for the abundance is solved numerically 



starting from large temperatures. However, this procedure poses a problem for Eq. (41). The step of 
the numerical solution is inversely proportional to A(T) and as long as A{T) is not suppressed by the 
Boltzmann factor included in Y, the step is too small and the numerical method fails. 

To avoid this problem, we use the fact that at large temperatures one can neglect the Q term in 



Eq. (41) and write the explicit solution for the linearised equation. The approximate solution in the 
case of large A is 

AY,{s) = A;^{s)C,{s). (47) 



One can use Eq. (47) to find the lowest temperature where AYi ~ 0.051^ and start solving numerically 
Eq. (41 ) from this temperature. In the general case it gives a reasonable step for the numerical solution 
Ss/s ~ 0.1, where s is the variable of integration. This method can, however, lead to some numerical 
problems if the masses of the two dark matter particles are very different. Let us call the light particle 
[ and the heavy particle ft. We have to start the numerical solution at a temperature T above the 
freeze-out temperature of the heaviest DM, 

Tfo^ « Ma/25. (48) 

At this temperature, 

Yt ^ M,-Mc 

— Riexp— ^ , (49) 

ifi J- iofi 

and the step in the numerical solution of the two component equations will be suppressed by a factor 
exp (—Mn — M[)/T{oii- This small step size is problematic when solving numerically the equation 
with the Runge-Kutta method. This occurs when M^/Mc > 2. In this case the equation for the 
heavy component must be solved independently assuming that the light component has reached its 
equilibrium density. If Mn/Mc < 2, the Runge-Kutta procedure can be used to successfully solve the 



thermal evolution equations (41) 



The abundances li and I2 will be modified by the interactions between the two dark matter 
sectors]^ Thus the new terms in Eq. (36) will simply add to the standard annihilation process with 



SM particles and will contribute to decrease the final abundance Yi. After X2 freezes-out, interactions 
of the type 22 — )■ 11 lead to an increase of Y2. When M^;^ <^ Mx2, the evolution of Y2 will be 
strongly infiuenced by the first sector since at its freeze-out temperature Yi is large. Following the 



same argument as above the new annihilation terms in Eq. (37) will contribute to a decrease in the 
final abundance Y2. Furthermore, the semi-annihilation process 12 — )• 10 which is always kinematically 
open means that xi acts as a catalyst for the transformation of X2 into SM particles. Thus the light 
component forces the heavy one to keep its equilibrium value, resulting in a significant decrease of 
the relic density of X2- When both DM particles have similar masses, the interplay between the two 
sectors is more complicated, in particular the role of the interactions of the type 20 — )• 11 will depend 
on the exact mass relation between the two DM particles. For example, this interaction can lead to 
an increase of the abundance of X2 if Yi is large enough for the reverse process to give the largest 
contribution. 



4.2 Numerical results 

The scalar model with a Z4 symmetry contains two dark sectors. In sector 1 the DM candidate is 
a complex singlet, S, the main contribution to crl^^^ comes from annihilation into Higgs pairs and is 
determined by the term AsijS'Pli^ip. Sector 2 is similar to the Inert Doublet Model (IDM). The DM 
candidate can be either the scalar or the pseudoscalar A^. Annihilation of DM into SM particles is 

*Note that Yi and Y2 correspond to the abundances of the particles with a given Z4 charge. The relative size of the 
masses of the DM particles depend on the choice of parameters in a given model. 
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usually dominated by gauge boson pair production processes, while annihilation into fermion pairs as 
well as co-annihilation processes can also contribute. Furthermore, for a DM mass at the electroweak 
scale, it was shown in [37| that annihilation into 3-body final states via a virtual W can be important 
below the W threshold. To avoid this complication we will consider a DM with a mass above masses of 
the W, Z, and h. Under this condition, the DM annihilation into SM particles in sector 2 is driven by 
SU{2) X ^7(1) gauge interactions and leads typically to a value of ^lh^ < 0.1, except for a DM heavier 
than about 500 GeV. The co-annihilation of H°, A^, H+ states increases . 

We will consider a benchmark point where both DM candidates S and have a mass near 350 
GeV. Other parameters are chosen so that semi-annihilation processes play an important role, while 
both components have comparable relic density and O/i^ = Clih'^ + Q,2h'^ = 0.1. In particular to have 
n2h^ ^ 0.05 requires the contribution of coannihilation processes - we therefore impose a small mass 



splitting Mfjo ~ Af^o, meaning that A5 will be small, see Eq. (22). Furthermore, a small value of A4 
also leads to a small mass splitting with the charged Higgs. Note that for small A5 and A4 the positivity 
condition on the potential, Eqs. ( 2|15 ) is easily satisfied. 



A2 


0.1 


Asi 


0.1 


A's 


0.1 


Ma 


341 GeV 


A3 


0.1 


As2 


0.3 


^5 


100 GeV 


Mh 


339 GeV 


A4 


0.01 


As'12 


0.13 


Mh 


125 GeV 


Ms 


350 GeV 


As 


0.1 


As2i 


0.13 











Table 3: Benchmark point for Z4. 



The results of the calculation of the relic density when including different terms in Eq. ( 36|37 ) is 
presented in Table [4j When only (co-)annihilation into SM particles are taken into account, the relic 
density of S is too high, while annihilation is much more efficient in Sector 2. Adding the interactions 
of the type of 1, 1 -f-)- 2,2 brings the value of Oi/i^ and $12^^ closer to each other. In our example 
the DM in sector 1 has weak interactions with SM particles, therefore Q.ih? is large when sector 2 is 
neglected. As a result of interactions with sector 2 particles the value for Q.ih? is significantly reduced. 
This effect was also observed for a DM model with a Z2 x Z2 symmetry [17| and was called the assisted 
freeze-out mechanism. Finally, when semi-annihilation processes are included, both Qih? and VL2h? 
decrease. Note that for this benchmark point, the cross section for DM elastic scattering on proton 



included terms of Eq(36 
^1100 ^2200 

4^00^ ^2200 and ^1122 
Ah 



37) 



lllb li2' 

0.24 0.041 

0.079 0.064 

0.050 0.051 



Table 4: Relic density of DM particles for the Z4 benchmark point. 

and neutron is 1.5(1.8) • 10^^ pb for the DM in sector 1 and 2, respectively. This is well below current 
exclusion limits of XENONIOO [34], as will be discussed at the end of this section. 

To examine more closely the interplay between the two DM sectors as well as the role of semi- 
annihilation in determining the DM abundance, we let Ms vary in the range 200-600 GeV and solve 
for the relic density by including new terms one by one. All other parameters are fixed to the value 
for the benchmark in Table IHl 

First we consider only the impact of annihilation processes, the results are displayed in Fig. [2] (left). 
When solving the evolution equation for the two DM independently, Q.ih? rises rapidly with Ms while 
Q.2h remains constant. Note that for the model under consideration we have that (T^^°° < erf °° hence 
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Figure 2: Effect of interactions between the two dark matter sectors (left) and of semi- annihilation 
(right) on r2i/i^(solid) and $12/1^ (dashed) as a function of M5. Left panel - Including only al^^^ and 
(T^^°° (black) as weh as cr^^^^, tr^^^^ (red). Right panel - Including only cr^^^° (green), only c7^^^° (red) 
as well as all semi-annihilations (blue), as a reference in black r2i/i^(solid) and (^2/1^ (dot) with only 
standard annihilation terms. Note that cr^^^" does not change Qih'^. 



flih"^ > 0,2^'^- The impact of o"^^^^ and o"^^^^ on the relic density depends on the relative masses 
of the scalar and doublet DM. The heavier DM candidate freezes out at a larger temperature than 
the lighter one, Tfo^ > Tfo/-. If the mass difference is large, this happens when the light one is at its 
equilibrium value. Thus the contribution of a^^^'' just adds to ct^^°°, leading to a decrease of the heavy 
DM abundance. Furthermore, after the light DM freezes out, interactions such as /t/t — )• H give an 
additional source of light DM, while the reverse reaction is suppressed by a Boltzmann factor. This 
effect can, however, be small when the heavy particles have a low density at this point. Thus the 
interactions between DM sectors 1 and 2 lead altogether to a decrease of the abundance of the heavy 
component and an increase of the light component. This is observed in the left panel of Fig. [2j In the 
region where Ms < Mfjo = 350 GeV, ^2^^ decreases while in the region Ms > Mjjo, ^2h'^ increases 
and vice- versa for Qih'^. Note that for large values of Ms, interactions with SM particles are weak so 
cr^^^ ^ crl^^^, leading to a large decrease in f^i/i^. When there is a small difference between the two 
DM particles, the freeze-out temperatures of both component are similar. The density of the heavy 
DM component has not yet decreased to its final value at the time the light component freezes out, 
thus the effect of /t/t -f-T- [[ interactions in increasing the abundance of the light component is more 
important. This is particularly noticeable when looking at the curve for 1^2/1^ in the region, where Ms 
is just above Mj^o = 350 GeV in Fig. [2] (left). This discussion, where we ignore the semi-annihilation 
terms, applies to models with As'12 = As2i = 0. In this case the Z4 symmetry is replaced with a Z2 x Z2 
symmetry. 

Next we consider the impact of semi-annihilation processes, ignoring the annihilation of pairs of 
particles from sector 1 to 2. The o"^^^" term does not affect Oi/i^ and works as a catalyst for 2 — )• SM 
transitions. This term has an effect only after the freeze-out of and its effect is stronger when 
Yi is large, see Eq. (37). Thus in the region Ms > M^o where the freeze-out of S occurs first (at 
a higher temperature), we find a roughly constant factor of suppression of Q.2h?'. As M5 decreases, 
its abundance Yi at the freeze-out of iT^ofi) will increase, thus the suppression of 1^2^^ is more 
important, see Fig. [2] Note that the suppression of Vl2h? for Ms > Mfjo is significantly larger than 
for the other semi-annihilation processes that we will discuss below. This is because the o"^^^^ term in 



12 




14- 12 10 S 6 4 2 



T (GeV) 

Figure 3: Temperature evolution of Yi (solid) and Y2 (dashed) with standard terms and the contribu- 
tion of o'l^'^^ for Ms = 260 GeV. Temperature evolution of Y2 with only standard terms (green/dashed), 
T is in GeV. 



Eq. 37 depends on Y^, which is large in this approximation. 

The second type of semi-annihilation process, 11 — )• 20 (or its reverse 20 — )• 11) leads to variations 
in the relic density of both DM components. If Ms > M^jo, the impact of crl^"^^ is very similar to the 
one discussed above for cr^^^^. For S, the heavy component, the overall annihilation cross section is 
increased, leading to a decrease in fix, illustrated by the blue curve in Fig.[2| For H^, the relic density 
increases because the process 11 — )• 20 is an additional source of sector 2 particles. This increase is 
even more important when both particles have similar masses - see the blue dashed curve in Fig. [2] 
when Ms = 260-350 GeV. To examine more closely the impact of the semi-annihilation in the region 
where the mass of both DM particles are similar, we compute the temperature evolution of Yi and Y2 
choosing M5 = 260 GeV. The result is displayed in Fig. |3j in particular comparing the evolution of 
Y2 with and without the contribution of cr^^^". For this choice of masses, the freeze-out of occurs 
when the abundance Yi = Yi is large, this means that the term 

1 1120 fv2 ^2 _ , 1 ^1120T72 f y2 



^:,-^n^-y^v1J=+-.-r;^^=^-i)... (50) 



in Eq. (37) forces Y2 to follow its equilibrium value. Thus I2 is further reduced by semi-annihilation at 
large temperatures. After the freeze-out of S, when Yi S> l^i, the same interaction leads to an increase 
of Y2 . Thus the overall effect is an increase in the abundance of class 2 particles as compared with the 
case where only standard interactions are considered. 

Finally, when Ms < 260 GeV, the cross section c^^^^, which consists of processes of the type 
SS — )• H^h is small because of a lack of phase space, thus Q.ih'^ is the same as when only standard 
annihilation terms were included. At the same time the reverse process, 20 — )• 11 drives the depletion 
of class 2 particles and ^2/1^ drops to very small values. Note that when M^jo > 2Ms we expect that 
the class 2 DM will decay into pairs of class 1 particles since they are allowed by the Z4 symmetry. 
However, in this example, the effect of o"^^^*^ and cr^^^^ terms already leads to very small values of 
02/1^ for low values of Ms, so that the decays are irrelevant. In summary, the combined effect of 
semi-annihilation processes is for this example close to the result of only including 

cjii20, see Fig.§ 

The result for Q.ih'^ and $72^^ including all annihilation and semi-annihilation processes is displayed 
in Fig. |4j The semi-annihilation mechanisms dominate for Ms < Mjjo while the assisted freeze-out 
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Figure 4: Left: f^i/i^ (solid), 1^2/1^ (dashed) and O/i^ (green) as a function of M5, the singlet DM mass. 
Right: Number of events expected in XENONIOO from S (solid) and (dashed) elastic scattering as 
a function of Ms- 



mechanism is the dominant effect when Ms > Mjjo . The total dark matter abundance is within about 
10% of the value preferred by WMAP measurements over the whole range of masses considered. While 
the features we have described here are generic, the relative importance of different annihilation and 
semi-annihilation processes is model dependent and depends on the size of the various cross sections 
within a specific model. 

Finally we compute the spin-independent cross section for S and scattering on xenon nuclei. 
As mentioned above at the benchmark point a^^ = 1.5(1.8) x 10^^ pb for the DM in sector 1 and 2 
respectively. We then compute the number of events that should be expected in XENONIOO [34j in 
the interval 8.4 keV < E < 44.6 keV after an exposure of 1171 kg-day. The number of events is directly 
proportional to the DM local density and we assume that the fraction of each DM component locally 
is the same as in the early universe, pi = pfij/Otot where p = 0.3. For S the cross section is largest for 
small masses, furthermore S contributes maximally to the DM density, hence the maximum predicted 
number of events, see Fig. H The cross section for scattering on nuclei is clearly independent of 
Ms, the variation of the number of events is simply due to the variation in the density of the second 
DM. 

5 CONCLUSIONS 

We have formulated scalar dark matter models with the minimal particle content in which dark matter 
stability is due to the discrete Z]sf symmetry with N > 2. Already the minimal models containing 
one extra scalar singlet and doublet possess non-trivial dark matter phenomenology. In particular, 
the annihilation processes with new topologies like XiXj — )• XkX, where Xi is one of the dark matter 
particles and X is any standard model particle, change the dark matter freeze-out process and must 
be taken into account when calculating the dark matter relic abundance. Furthermore, in models 
with two dark matter candidates, annihilation processes involving only particles of two different dark 
matter sectors also impact the relic abundance of both dark matter particles. We have performed an 
example study of semi-annihilations in two scalar dark matter models based on Z3 and Z4 symmetries. 
We implemented those models for micrOMEGAs and studied the impact of semi-annihilations and 
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of the interactions between the dark sectors on the generation of dark matter rehc abundance at the 
early Universe and the predictions for dark matter direct detection relevant for the presently running 
XENONIOO experiment. We conclude that in this type of models both semi-annihilations and dark 
sector interactions may significantly affect the dark matter phenomenology compared to the well studied 
Z2 models, and, therefore, must be taken into account in precise numerical analyses of dark matter 
properties. 
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